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ABSTRACT 

The GRB 991216 and its relevant data acquired from the BATSE experiment and RXTE 
and Chandra satellites are used as a prototypical case to test the theory linking the origin of 
gamma ray bursts (GRBs) to the process of vacuum polarization occurring during the formation 
phase of a black hole endowed with electromagnetic structure (EMBH). The relative space-time 
transformation paradigm (RSTT paradigm) is presented. It relates the observed signals of GRBs 
to their past light cones, defining the events on the worldline of the source essential for the 
interpretation of the data. Since GRBs present regimes with unprecedently large Lorentz 7 
factor, also sharply varying with time, particular attention is given to the constitutive equations 
relating the four time variables: the comoving time, the laboratory time, the arrival time at the 
detector, duly corrected by the cosmological effects. This paradigm is at the very foundation of 
any possible interpretation of the data of GRBs. 

Subject headings: black holes, gamma ray bursts, supernovae 



In recent years, a large variety of very accurate 
experimental data, ranging from 7 rays all the way 
to the radio band, has been obtained for the after- 
glows of GRBs, following their first discovery by 
the BeppoSAX satellite (see e.g. Costa 2001, and 
references therein). 

In the theoretical models of GRBs there are cur- 
rently three topics under debate: 
1) The "internal Shock Model", introduced by 
Rees & Meszaros (1994), many aspects of which 
have been developed (Paczyhski & Xu 1994; Sari 
& Piran 1997; Fenimore 1999; Fenimore et al. 
1999). The underlying assumption of this model 
is that all the variations of GRBs in the range 
At ~ 1 ms up to the overall duration T of the or- 
der of 50 s are determined by the "inner engine". 
The difficulties of explaining the long time scale 

1 ruffini@icra.it 



bursts by a single explosive model has evolved into 
a class of models assuming an "inner engine" with 
a prolongued activity (see e.g. Piran 2001, and 
references therein). 

2) The "external shock model", also introduced by 
Meszaros & Rees (1993), is less popular today. It 
relates the GRBs' light curves and time variations 
to interactions of a single thin blast wave with 
clouds in the external medium. There is the dis- 
tinct possibility, within this model, that "GRBs' 
light curves are tomographic images of the den- 
sity distribution of the medium surrounding the 
sources of GRBs" (Dermer & Mitman 1999) (see 
also Dermer et al. 1999; Dermer 2000, and ref- 
erences therein). In this case, the structure of the 
burst does not come directly from the "inner en- 
gine" . 

3) In order to decrease the energy requirements of 
GRBs, the effect of beaming has been advocated 
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(see e.g. Mao & Yi 1994; Davies et al. 1994). The 
possibility of inferring its existence from changes 
in the power-law index of the afterglow is gen- 
erally considered attractive (see e.g. Meszaros & 
Rees 1997a; Rhoads 1997; Meszaros et al. 1998 
Panaitescu et al. 1998; Dermer & Chiang 1999 
Sari at al. 1999; Panaitescu & Meszaros 1999 
Rhoads 1999; Halpern et al. 2000; Gou et al. 
2001). 

For the astrophysical nature of the system orig- 
inating the GRB, a binary system of merging neu- 
tron stars has been proposed (see e.g. Eichler et 
al. 1989; Narayan ct al. 1992; Meszaros & Rees 
1992a, b). Problems occur: a) in the general en- 
ergetics which cannot exceed ~ 3 x 10 52 ergs, b) 
to explain the longer bursts (see Salmonson et al. 
2001; Wilson et al. 1996), and c) in the observed 
location of the GRBs' sources in star forming re- 
gions (sec Bloom et al. 2000). Alternatively novel 
classes of astrophysical systems have been postu- 
lated, including black hole - white dwarf (Fryer et 
al. 1999) and black hole - neutron star binaries 
(Paczyhski 1991; Meszaros & Rees 1997b), as 
well as Hypernovae (see Paczyhski 1998), failed 
supernovae or collapsars (see Woosley 1993; Mac- 
Fadyen & Woosley 1999) and supranovae (see Vi- 
ctri & Stella 1998, 1999). 

We take a somewhat intermediate approach by 
studying the GRB emitted by the process of vac- 
uum polarization around a black hole endowed 
with electromagnetic structure: the EMBH model. 
Such a model has the advantage that all the ba- 
sic intermediate theoretical background, starting 
from the process of gravitational collapse itself, 
have been developed. The model can therefore 
make precise predictions which can be compared 
with the observations. 

We consider a GRB "prototypical" case, in or- 
der to create a new interpretative paradigm, to be 
then applied to the observations of other GRBs. 
Since some of the best data available, from BATSE 
(BATSE Rapid Burst Response 1999), RXTE 
(Corbet & Smith 2000) as well as the remark- 
able accuracy of the Chandra (Piro et al. 2000) 
satellite are available for the GRB 991216, we use 
it as the prototype. In addition: a) it is one of the 
strongest observed GRBs. b) it radiates mainly in 
X- and 7-rays and less than 3% is emitted in op- 
tical and radio band, and c) a precise value of the 
slope of the energy emission during the afterglow 



as a function of time, n = —1.64 (Takeshima et 
al. 1999) and n = -1.616 ± 0.067 (Halpern et al. 
2000), has been obtained. 

The EMBH model relates the origin of the en- 
ergy of GRBs to the extractable electromagnetic 
energy of an EMBH (Christodoulou & Ruffini 
1971) via the vacuum polarization process oc- 
curring during the gravitational collapse leading 
to the formation of an EMBH (Damour & Ruffini 
1975). The first step in this theory is the definition 
of the dyadosphcre (Ruffini 1998; Preparata et 
al. 1998), an extended region outside the EMBH 
horizon formed of an optically thick plasma of 
electron-positron pairs and radiation whose en- 
ergy Edya is related to the mass fi = M/M© and 

electromagnetic parameter £ = Q/ (MVG^j of the 
EMBH by the relation: 
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where r+ = 1.47 x 10 5 Ai(l + ^1 - £ 2 ) is the 
horizon radius and rd s = 1.12 x lO 8 ^/^ is the 
dyadosphere radius and, as usual, M and Q are 
the mass-energy and charge of the EMBH and G 
is the Newton constant of gravity. 

The evolution of this pair-electromagnetic 
plasma leads to the formation of a sharp pulse (the 
PEM pulse) that very rapidly reaches a Lorentz 
gamma factor of 10 2 and higher. The subsequent 
interaction of this pulse with the baryonic matter 
of the remnant, left over from the gravitational 
collapse of the protostar, and with the interstel- 
lar medium (ISM) leads to the different eras of 
the GRBs. It is useful to parametrize the bary- 
onic mass Mb of the remnant by introducing the 
dimensionless parameter B: 



M B c 2 = BE dy 



(2) 



The confrontation of the theoretical model with 
the observational data allows an estimate for the 
values of the EMBH parameters. It also allows us 
to probe the density of the baryonic material in 
the remnant, in the ISM as well as in the stellar 
distribution within a few parsecs of the EMBH 
(see Ruffini et al. 2001b,c). 

The first step in this process is the establish- 
ment of the first set of constitutive equations re- 
lating: 
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a) The comoving time of the pulse (r): the evolu- 
tion of the thermodynamical quantities (density, 
temperature) are computed using this time. 

b) The laboratory time (t) defined by an inertial 
reference frame in which the EMBH is at rest. 

c) The arrival time (t a ): the laboratory time at 
which light signals from the source reach a distant 
observer at rest in the laboratory frame. The zero 
of the arrival time has been chosen to coincide with 
the arrival of the light signals from the moment of 
formation of the EMBH 

d) The arrival time at the detector (if ): this is the 
arrival time taking into account the cosmological 
redshift of the GRB source. We have 



f" 



ta (1 + Z) , 



(3) 



where z is the cosmological redshift of the GRB 
source (Ruffini ct al. 2001f; Bianco ct al. 2001c). 
In the case of GRB 991216 we have z ~ 1.00. 

The mutual relations of these four times with 
the radial coordinate in the laboratory frame is 
the subject of this letter. We first give emphasis 
to a basic feature of the arrival time determina- 
tion. For signals emitted by a pulse moving with 
velocity v in the laboratory frame, we have: 



At a = [t + At 



to 



^)=At- r -, 
c I c 



(4) 

where At a (At) is the time interval in arrival (lab- 
oratory) time, i?o is the distance of the observer 
from the EMBH, to is the laboratory time cor- 
responding to the gravitational collapse, and r 
is the radius of the expanding pulse at the time 

t = t + At. 

For simplicity we abbreviate the interval nota- 
tion At a (At) by t a (t). Eq.(4) can then be rewrit- 
ten as: 



r J v(t')dt' 

t a =t--=t- J0 

C C 
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where the dyadosphere radius Tds is the value of 
r (t = 0). We consider only the photons emitted 
along the line of sight, since the spreading due to 
the angular dependence and to the thickness of the 
pulse is negligible (see Bianco et al. 2001b, for 
details). The solution of Eq.(5) neglecting r^s- 
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is in general highly nonlinear. 

If and only if v is constant and v 
can be rewritten, neglecting rd s , as: 



c, Eq.(5) 
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It is clear that the knowledge of t a , which is in- 
deed essential for any physical interpretation of 
GRB data, depends on a definite integral whose 
integration limits extend from the gravitational 
collapse to the time t relevant for the observa- 
tions, see Eq.(5). Such an integral is not generally 
expressible as a simple linear relation or even by 
any explicit analytic relation since we are dealing 
with processes with variable Lorentz gamma fac- 
tors of unprecedented magnitude and time vari- 
ability. Most studies have adopted an approxi- 
mation of the kind given in Eq.(7) (see e.g. Fen- 
imore et al. 1996). We instead use Eq.(5). The 
adoption of Eq.(7) misses a crucial feature of the 
GRB process and it leads to a subversion of the 
space-time relations in GRBs, with a wide range of 
consequences: all theoretical computations on the 
power-law indexes of the afterglow are affected. 
Specific illustrative examples, pointing out these 
differences, are shown in the next paragraphs (see 
Ruffini et al. 2001f, for details). 

The book-keeping of the four different times 
and corresponding space variables must be done 
carefully in order to keep the correct causal rela- 
tion in the time sequence of the events involved. 
This will have also important consequences in the 
supernovae-GRBs correlation (see Ruffini ct al. 
2001c) 

The second set of constitutive equations are 
the full non-linear relativistic hydrodynamic equa- 
tions of energy and momentum conservation, to 
be solved together with the rate equation for the 
plasma. The computations carried out scmi- 
analytically in Rome have been validated by the 
full numerical computations performed using Wil- 
son's codes at Livermore (see Ruffini et al. 1999, 
2000, 2001f). 

We have integrated both sets of constitutive 
field equations given in Ruffini et al. (1999, 2000); 
Bianco et al. (2001a) , for the source GRB 991216. 
Correspondingly, we have obtained the parame- 
ter values presented in Ruffini et al. (2001b): 
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E dya ~ 10 53 ergs and B = 4x 10~ 3 . These values 
correspond to any of the following pairs of values 
for the EMBH mass and charge to mass ratio (/x, £) 
= (22.3, 0.1); (10.0,0.15); (5.5,0.2). 

Crucial to any GRB data interpretation is the 
relation of the Lorentz gamma factor to the radial 
coordinate of the source in the laboratory frame 
and the corresponding values of the above four 
time parameters. In Fig. 1 the gamma factors 
for the different eras are given as a function of 
the radial coordinate of the source in the labora- 
tory frame. Correspondingly we present in Fig. 2 
the relation between the laboratory time and the 
detector arrival time for the source GRB 991216. 
The highly nonlinear behaviour is obvious, and the 
different results obtained from the use of Eqs.(6,7) 
are clearly visible. Details are given in Ruffini et 
al. (2001f). 

In Tab. 1, for each successive "era" and for one 
very significant event, we give the initial- final val- 
ues of the gamma Lorentz factor, the four time 
parameters mentioned above, as well as the radial 
coordinates in the laboratory frame. 

We then have: 

1) Era I. The pair-electromagnetic plasma, ini- 
tially at 7 = 1, expands away from the EMBH 
horizon and from the dyadosphere as a pulse, the 
PEM pulse. In the comoving frame the thickness 
of the pulse increases during the expansion, but 
the Lorentz contraction in the laboratory frame 
exactly balances this expansion so that, in the 
laboratory frame, a constant thickness approxima- 
tion can be adopted for the burst (Ruffini et al. 
1999). The expansion of the PEM pulse occurs in 
a region of very low baryonic contamination with 
density p B < 10~Vcto 3 (Ruffini 2001). The fi- 
nal Lorentz factor and space time parameters are 
given for point 2 in Tab. 1. 

2) Era II. While the PEM pulse is still opti- 
cally thick it reaches the remnants left over by 
the gravitational collapse of the progenitor star. 
The cngulfmcnt of this baryonic material induces 
by conservation of energy and momentum a dras- 
tic reduction in the 7- factor (Ruffini et al. 2000). 
The amount of baryonic matter in the remnant has 
been fixed by the determination of the parameter 
B in the fitting of the afterglow data (see Ruffini 
et al. 2001b). Since these data contain impor- 
tant direct information on the progenitor star, we 



report in Tab. 2 some specific values of the pa- 
rameters corresponding to selected values of the 
EMBH masses: they include the radius and thick- 
ness of the remnant as well as the density of bary- 
onic matter. The results are largely independent 
of the thickness A of the remnants (Ruffini et al. 
2000); they depend crucially only on the value B. 
The final Lorentz factor and space-time parame- 
ters are given for point 3 in Tab. 1. 

3) Era III. A new pulse is formed composed of 
electron-positron pairs and baryons and electrons 
of the remnant material, a PEMB pulse. Since 
the opacity from Thomson scattering consequently 
increases, the process of self-acceleration of the 
burst continues to even larger values of the Lorentz 
gamma factor, which may reach values up to 10 3 - 
10 4 in some sources (Ruffini et al. 2000). In the 
present case of GRB 991216, the maximum value 
reached is 239.6. It is remarkable that the con- 
stant thickness approximation for the pulse is still 
valid (Ruffini et al. 1999, 2000). This era ends at 
point 4 as the condition of transparency is reached. 
At that point, what we define here as the proper 
gamma-ray-burst (P-GRB) is emitted (Ruffini et 
al. 1999, 2000; Bianco et al. 2001a). This defini- 
tion is assumed in order to distinguish the overall 
GRB phenomena from the specific emission oc- 
curring as the moment of transparency is reached. 
The shape, the intensity and the duration of this 
P-GRB, also called in the literature the elemen- 
tary spike (Bianco et al. 2001a), strongly de- 
pend on the baryonic matter content (Bianco et al. 
2001b; Ruffini et al. 2001f). The final Lorentz fac- 
tor and space-time parameters are given for point 
4 in Tab. 1. 

These first three eras have no counterparts in 
earlier models, since no detailed description of the 
early phases of GRB has been attempted. 

4) Era IV. The accelerated baryonic mat- 
ter expands still as a constant thickness pulse 
(ABM pulse) at ultra-relativistic velocities and 
engulfs baryons and electrons from the interstel- 
lar medium (ISM), which is assumed to have a 
constant number density rii sm of 1 proton/cm 3 . 
We have assumed that all internal energy made 
available by the relativistic conservation of energy 
and momentum is radiated away in the afterglow, 
mainly in 7 (~ 90%) and X-rays (<~ 10%), and a 
few percent in the optical and radio emission (see 
Halpern et al. 2000). We have used the "fully 
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radiative case" condition (see e.g. Piran 1999), 
considering only the leading contribution of the 
head-on flux. We have neglected the spreading 
due to off-axis emission, considered e.g. in Feni- 
more et al. (1999) and have given physical rea- 
sons for neglecting such contributions (Ruffini et 
al. 2001f). At any specific time, the total flux, 
and consequently the bolometric luminosity of the 
afterglow, is fixed by the above requirements. The 
detailed spectral distribution depends on the dom- 
inant radiative process of the internal energy. A 
variety of such process have been considered in the 
literature and their results within the context of 
our model are given in Ruffini et al. (2001f). In 
Fig. 1 we show how, during this era, the Lorentz 
gamma factor first coasts to a constant value and 
then rapidly decreases, going from 7 = 239.6 to 
7 = 2.7. Most important is that the point P where 
7 ~ 160.2 corresponds to the peak of the afterglow 
(see Ruffini et al. 2001b). Beyond this point P the 
slope of the afterglow flux, as a function of arrival 
time, approaches the power law index n = —1.6 in 
perfect agreement with the observations of RXTE 
and Chandra (see Ruffini et al. 2001b). The fi- 
nal Lorentz factor and space-time parameters are 
given for point 5 in Tab. 1. It is important to em- 
phasize that this power law index results from the 
combination of three critical assumptions: a) the 
emission occurring in a "fully radiative" regime, 
b) the condition of spherical simmetry, and c) the 
constancy of the ISM density. Earlier results rele- 
vant to this treatment can be found in Sari (1997) 
and in Dermer et al. (1999) (see, for comparison 
and contrast, Ruffini et al. 2001f). 

5) Era V. This is the transition to the relativis- 
tic and nonrelativistic regimes. This era is more 
complex. It contains two successive sub-eras, one 
with a power-law index of the energy emitted in 
the afterglow as a function of the detector arrival 
time n = —1.36, corresponding to a still relativis- 
tic era (1.1 < 7 < 2.7), and a final one approach- 
ing the pure newtonian regime, with n = —1.45 
and 7 < 1.1. Comparison and contrast with exist- 
ing slopes in the literature (see e.g. Vietri 1997) 
are presented in Ruffini et al. (2001f). No data of 
GRB 991216 are available for checking the thorct- 
ical predictions of this last era. 

In conclusion we see from Tab. 1 and Fig. 1 the 
remarkable and perfectly reasonable results that a 
motion of the pulse corresponding to a displace- 



ment of 9.692 x 10 13 cm will correspond to an ar- 
rival time interval of 1.360 x 10 _1 sec, leading to 
what has been called apparent superluminal be- 
havior. Similarly, on a larger scale, a displacement 
of the pulse by 2.958 x 10 17 cm will correspond to 
an increment of 2.439 x 10 5 sec in arrival time, 
leading again to apparently superluminal behav- 
ior. 

From the above results we are ready to express 
the relative space-time transformation (RSTT) 
paradigm: the necessary condition in order to in- 
terpret the GRB data, given in terms of the arrival 
time at the detector, is the knowledge of the en- 
tire worldline of the source from the gravitational 
collapse. In order to meet this condition, given 
a proper theoretical description and the correct 
constitutive equations, it is sufficient to know the 
energy of the dyadosphere and the mass of the 
remnant of the progenitor star (see Ruffini et al. 
2001b). The application of this RSTT paradigm 
will have important consequences for the interpre- 
tation of the burst structure (IBS), leading to a 
new paradigm (see Ruffini et al. 2001b), as well 
as for the GRB-supernova correlation (Ruffini et 
al. 2001c). 

We thank three anonymous referees for their 
remarks, which have improved the presentation of 
this letter. 
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Fig. 1. — The theoretically computed Lorentz gamma factor for the parameter values Ed ya = 9.57 x 10 52 erg, 
B = 4x 10~ 3 is given as a function of the radial coordinate in the laboratory frame. The corresponding values 
in the comoving time, laboratory time and arrival time are given in Tab. 1. The different eras, indicated by 
roman numerals, are illustrated in the text, while the points 1,2,3,4,5 mark the beginning and end of each 
of these eras. The point P marks the maximum of the afterglow flux (see Ruffini ct al. 2001b). At point 4 
the transparency condition is reached. 
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Fig. 2. — Relation between the arrival time (tf) measured at the detector, and the laboratory time (t) 
measured by an observer at rest in an inertial frame also at rest with the GRB source. The solid curve 
is computed using the exact formula given in Eq.(5), incorporating as well the cosmological effects given 
in Eq.(3), while the dashed-dotted curve is the corresponding line obtained using the approximate formula 
given in Eq.(7), computed with the time varying 7 given in Fig. 1. This is the approximation used by all 
recent works on GRBs. The difference is very conspicous, if one takes into account that the diagram is in 
logaritmic scale, and has basic consequences on the astrophysical scenario (see e.g. Ruffini et al. 2001b, c). 
We distinguish four different phases. Phase A: There is a linear relation between t and (dashed line). 
Phase B: There is an "effective" power-law relation between t and t% (dotted line). Phase C: No single 
analytic formula holds, the relation between t and can only be analysed patchwise and has to be directly 
computed by the integration of the complete equations of energy and momentum conservation. Phase D: 
As the Lorentz factor approaches 7=1, the relation between t and if asymptotically goes to t = t% (light 
gray line). Details given in Ruffini et al. (2001f). 
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Table 1 

lorentz factors for selected events and their space-time coordinates. 



Point 


r (cm) 


r(s) 


t(s) 


ta (S) 


4 (s) 


7 


1 


1.610 x 10 8 


0.0 


0.0 


0.0 


0.0 


1.00 


2 


7.659 x 10 9 


1.985 x 10~ 2 


2.580 x lO" 1 


1.846 x 10~ 3 


3.692 x 10~ 3 


48.38 


3 


9.153 x 10 9 


2.292 x 10~ 2 


3.089 x 10- 1 


2.780 x 10~ 3 


5.559 x 10~ 3 


11.38 


4 


9.692 x 10 13 


14.23 


3.295 x 10 3 


6.805 x 10~ 2 


1.361 x lO" 1 


239.6 


P 


4.863 x 10 16 


7.784 x 10 3 


1.653 x 10 6 


11.86 


23.72 


160.2 


5 


2.958 x 10 17 


1.082 x 10 6 


9.989 x 10 6 


1.2195 x 10 5 


2.439 x 10 5 


2.7 



Table 2 

baryonic matter of the remnant. 



M a 




r ds b 


r b 
' shell 


A b 

^shell 


P c 


22.3 


0.10 


1.67 x 10 8 


8.36 x 10 9 


1.67 x 10 9 


0.30 


10.0 


0.15 


1.37 x 10 8 


6.86 x 10° 


1.37 x 10 9 


0.55 


5.5 


0.20 


1.17 x 10 8 


5.8 x 10 9 


1.17 x 10 9 


0.90 



a in solar mass 

b in the laboratory frame in cm 
c in g.cm~ 3 
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